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Using the full-potential linearized augmented plane-wave method within the density functional 
theory, we compute all contributions to the scattering independent part of the thermoelectric con- 
ductivity tensor axy, namely the intrinsic contribution a^^y and the side-jump contribution a^iy. For 
the ferromagnetic materials bcc Fe, hep Co, fee Ni and Llo ordered alloys FePd and FePt, our in- 
vestigations of the energy and temperature dependence of the intrinsic and side-jump contributions 
show that they are both of equal importance. Overall, our calculations are able to correctly repro- 
duce the order of magnitude and sign of the experimentally measured signal, suggesting that the 
scattering independent part plays an important role in the anomalous Nernst effect of ferromagnets. 

PACS numbers: 



When a thermal gradient is applied to a ferromagnetic 
material that is placed in a magnetic field, the anomalous 
Nernst effect (ANE) shows itself as an electric field E 
that emerges in a direction perpendicular to the applied 
thermal gradient — VT and perpendicular to the sample's 
magnetization M.^ Prominently, the Nernst signal can 
be used as a probing tool for the vortex phase in type 
II superconductors^''^ and it has been also discussed in 
spinel ferromagnets"* and on a surface of a topological 
insulator.^ However, despite increasing interest in this 
phenomenon in the past years, ^ no attempt to predict 
from first principles the values of the ANE in metallic 
ferromagnets, such as bcc Fe, has been made so far. 

It is convenient to quantify the ANE in terms of 
the thermoelectric conductivity tensor a, also called the 
Peltier or Nernst conductivity tensor. In linear response 
theory, the expression for the charge current j in the pres- 
ence of a weak electric field and a small thermal gradient 
reads'' 

j = (T • E - a • VT, (1) 

where the electric conductivity tensor is denoted by a. 
The matrix elements of a and a are related via the gen- 
eralized Mott formula'^"^*^ 

axy^-- j de — a^y (2) 

where e = — |e| is the electronic charge, e the energy and 
^ the chemical potential of the electrons which appears in 
the Fermi distribution function /(e,/i,T). In ferromag- 
netic materials, we can replace the off-diagonal conduc- 
tivity axy with the anomalous Hall conductivity cr^^^,^^ 
which is usually decomposed into an intrinsic Berry cur- 
vature driven contribution crj^^ and two extrinsic, stem- 
ming from disorder, contributions. Of the latter two, the 
first is the so-called side-jump contribution cr^^y, which is 
caused by the scattering of electrons off impurities but 
which paradoxically does not depend on their concen- 
tration n.i}'^ The second disorder-driven contribution is 
the skew-scattering (t|'^, which is inversely proportional 



to the impurity concentration, i.e., cr^^ cx l/rii}^'^^ The 
Mott relation, Eq. (2), implies that the thermoelectric 
conductivity tensor can be decomposed in an analogous 
way: 

axy = aL'y + ot% + otfy (3) 

The difficulties in understanding the origins of the a- 
nomalous Hall effect (AHE) have greatly impeded the 
progress in the field of the ANE as well. For exam- 
ple, on the side of qualitative theory of solids, we are 
aware of only a single estimate for the value of a^^y in 
the cuprate CuCr2Se4_2;Bra;.^^ However, as we have re- 
cently shown, all scattering independent contributions to 
the AHE, that is, a^y and a^jy, can be calculated from 
first principles on an equal footing from the knowledge of 
the electronic structure of the pristine crystal alone, al- 
lowing for a quantitative comparison between theory and 
measurements performed on samples in the good metal 
regime, for which skew-scattering is suppressed.*® In the 
present work, we extend this methodology to the ANE. 
We calculate all the scattering-independent contributions 
to the thermoelectric conductivity tensor, a^^y and a^jy, 
in bcc Fe, hep Co, fee Ni and LIq ordered alloys FePd 
and FePt. By comparison to experimental data, we show 
that a^^y and a^^y provide the correct order of magnitude 
and sign of the anomalous Nernst signal in transition- 
metal ferromagnets. We also make predictions concern- 
ing the temperature dependence of the scattering-inde- 
pendent ANE. 

Our approach is based on electronic structure calcula- 
tions performed within the full-potential linearized aug- 
mented plane- wave method as implemented in the Jiilich 
density functional theory code FLEUR.*^ The matrix el- 
ements of the multi-band Bloch Hamiltonian 77(k) in 
the basis of maximally-localized Wannier functions^"'^* 
have been computed using the Wannier interpolation 
technique^^ and inserted into the equations for a^^y and 
CT^jy obtained within the Kubo-Stfeda formalism, ^'^ as- 
suming short-range disorder in the system.^** The knowl- 
edge of any free parameters is not required in this scheme. 
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FIG. 1: (Color online) ANE at T = 300 K in Fe, Ni, Co, FePd and FePt. The bar diagrams in the first row of (a)-(e) show 
the values of a^^y and a^y at different energy positions. The true Fermi energy level in each material lies at the position of 
the middle vertical line. The line diagrams in the lower row of (a)-(e) depict the energy dependence of a]^y and a%^y. A solid 
pattern for bars and lines stands for the [001] magnetization direction, hatched bars and dotted lines for the [110] magnetization 
direction in Co, FePd and FePt. For Ni, also the LDA+C/ calculations with U = 3.9 eV and J — 2.6 eV are shown in (c). A 
sketch of the experimental setup is shown in (f). 



For the evaluation of the integral in Eq. (2), we adopted 
an energy grid that was denser at low temperatures, since 
the energy derivative of the Fermi function becomes a S- 
distribution in this limit. Near room temperature, we 
found that an energy spacing of Ae 5 meV offered the 
best trade-off between accuracy and computational cost. 

If follows from an argument by Berger^^ that large val- 
ues of axy could not arise if electrons of different energies 
experienced an AHE of the same magnitude and sign, 
since the transverse velocities of electrons diffusing down 
the applied temperature gradient would then cancel with 
those of the less energetic electrons diffusing up the tem- 
perature gradient and the net transverse current would 
be zero. The same can also be deduced from the Mott 
formula Eq. (2): Since the product (9^/) ■ {e — fi) is an 
antisymmetric function with respect to the Fermi energy 
level ep = A*, the ANE would vanish if the AHE conduc- 
tivity cTxy^^ a symmetric function around sp, i.e., if 



it was equal for a pair of energy values with the same 
distance to e^. 

In Fig. l(a)-(e), the values of the intrinsic and side- 
jump contribution to the AHE and the values of the ANE 
at room temperature in Fe, Ni, Co, FePd and FePt cal- 
culated from Eq. (2) as a function of the Fermi energy 
are shown. In bcc Fe, we observe that a}^y is nearly sym- 
metric around the Fermi level of 9.71 eV, resulting 
in a rather small value of a!f„. If we shift the Fermi en- 
ergy artificially by -1-0.08 eV, the course of a^^y gets more 
asymmetric with respect to the new Fermi energy and 
the value of a]^y increases significantly by a factor of 6. 
As a consequence, the total value a]^y + a^y changes sign 
and becomes positive. If the Fermi energy level is low- 
ered to 9.6 eV instead, all contributions become negative. 
In general, we find that the thermoelectric conductivity 
in Fe, Co and Ni appears to be highly sensitive to the 
position of the Fermi level, suggesting that the ANE in 
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ferromagnets can be easily tuned by e.g. suitable dop- 
ing. On the other hand, for the compounds FePd and 
FePt the Fermi energy dependence of axy is much less 
pronounced. 

By now it is established that the side-jump contribu- 
tion to the AHE is important in FePd whereas the in- 
trinsic AHE is dominant in FePt.^^'^^ As follows from 
Fig. 1(d)- (e), this statement also applies to the ANE in 
these materials, i.e., a^jylep) is as large as a^^y{eF) in 
FePd, but only half of this value in FePt. This crossover 
behavior is caused by the different spin-orbit interaction 
strength of Pd and Pt atoms. In Fe and Co, the mag- 
nitude of a^y{eF) is greater than that of a^^y{eF), albeit 
(Jl^yisp) being smaller than cr^^yiep) in both materials. 

In analogy to the AHE,"'*'^ one might suspect that the 
ANE is highly anisotropic with respect to the direction 
of the magnetization in the crystal for hep Co, and LIq 
ordered FePd and FePt alloys, due to their uniaxial crys- 
tal structure. Indeed, at the Fermi energy, the side-jump 
contribution a^^y switches its sign in Co and is strongly 
reduced in FePd and FePt as the direction of the magneti- 
zation is changed from [001] to [110] direction. However, 
the anisotropy of the intrinsic contribution al^y is not that 
pronounced. Such a different dependence of a'^y and a^Jy 
on the magnetization direction may be attributed to the 
different distribution of (T^^fy and a^^y in the Brillouin zone 
of these materials.^'' 

The ANE in fee Ni presents an exceptional case, since 
the theoretical value of the thermoelectric conductivity 
in this material is much larger than in other considered 
compounds, see Fig. 1(c). In Ni, the intrinsic anomalous 
Hall conductivity is sharply peaked near the Fermi en- 
ergy, and the respective value for the ANE depends on 
which side of the peak it is evaluated. There is profound 
evidence, however, that the large value of a}^y is mainly 
an artifact of the local density approximation (LDA) or 
the generalized gradient approximation (GGA), because 
correlation effects among the Bd electrons in Ni become 
of crucial importance for its properties. ^^■'^^'^^ For this 
reason, we have adopted the same approach as in our 
previous work and took the correlation effects into ac- 
count within the GGA+U scheme. '^'^ For the intra-atomic 
Coulomb repulsion and exchange parameters U and J 



TABLE I: Comparison of Q'^'^'^' with experimental data for 
Qi [10~'^^V/(K-G)] near room temperature. 
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we chose the values of 3.9 and 2.6 eV, respectively. This 
choice of parameters has been found to greatly improve 
the agreement of the calculated AHE in Ni, while keep- 
ing the magnetic moment close to the value observed in 
experiment. Fig. 1(c) reveals that correlations have a sig- 
nificant effect on the energy dependence of o'xy^^ axy 
in Ni. We observe that the peak in the intrinsic contri- 
bution to the AHE flattens out, whereas the side-jump 
contribution remains mostly unaffected upon including 
the U. This can again be understood from the different 
behavior of the two effects at the Fermi surface. Upon 
including U, the intrinsic contribution a]^y and the total 
value for a^^y + a^y changes its sign and the magnitude 
of the thermoelectric conductivity is greatly reduced. 

For comparison with experiment, we consider the 
situation in which a temperature gradient — VT — 
-~{dxT,dyT,dzT) II a; is applied to an electrically iso- 
lated sample perpendicular to the magnetization M || z, 
see Fig. 1(f). As a function of the magnetic induction 
B = H -I- 47rM, the Nernst effect obeys a law of the type 
Ey/dxT — H^Qo+^^T^MzQi, where Qq and Qi are the or- 
dinary and anomalous Nernst coefficients, respectively."^^ 
However, in ferromagnetic materials, the ordinary Nernst 
coefficient is very small, Qo Qi- The remaining co- 
efficient Qi is generally measured in a zero-current con- 
figuration, jx = jy = 0, with the boundary condition 
dyT = 0.^^ In this scenario, we obtain from Eq. (1) 

AttM^Qi = Pxxi^xy - SUxy) (4) 
= Pxxia'xy - S(Jxy) + Rxxi^xy " Sd^'^y), 

where Pxx = 1 /(^xx is the resistivity of the sample and the 
so-called Seebeck-coefficient is defined by S" = Ex/dxT = 
dxx/'^xx- The last line of Eq. (4) may be interpreted 
as AnMz{Qf + Q^i)^ where the intrinsic contribution to 
the anomalous Nernst coefficient is denoted by Q'f and 
the side-jump contribution is denoted by Q^^. While the 
Seebeck-coefficient S describes the conversion of a ther- 
mal current into a longitudinal electrical current, the 
Nernst-coefficient Qi is a measure of the corresponding 
transverse effect. Even though the value of the tempera- 
ture gradient —dxT and the strength of the magnetic in- 
duction Bz do not appear in Eq. (4) , they seem to have a 
great influence on the ANE experimentally.^^ Especially 
should be large enough that Mz is at the saturation 
magnetization value. However, values for Qi, pxx^ ctxy, 
S, axy and Mj, were not measured simultaneously in most 
experiments. We therefore gathered values for the resis- 
tivity Pxx and Seebeck coefficient S from various sources 
and computed the scattering-independent contribution 
Q'j^^^^ = Qf + Q^^ to the anomalous Nernst coefficient 
following Eq. (4) (Details of this procedure can be found 
in the supplement ).^^ 

A comparison of our calculated values with experimen- 
tal data is presented in Tab. (I). The experimental values 
show a considerable spread, which refiects the fact that 
the ANE is found to depend sensitively on experimental 
details and material-specific parameters. Nevertheless, it 
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FIG. 2: (Color online) Temperature dependence of the 
intrinsic- and side-jump contribution to the ANE. The val- 
ues of a'^y at T = are computed from Eq. (5). 



can clearly be seen that the inclusion of the side-jump 
contribution is crucial for Fe and FePt and lets theory 
and experiment match very well: Q^^^''^ is about 112% 
of the smaller experimental value in Fe and 105% of the 
experimental value in FePt. For Co, the side-jump con- 
tribution to Q^i^''^ is relatively small and only around 
10% for both magnetization directions, but it still brings 
the theoretical value for the anomalous Nernst coefficient 
closer to experiment. The lack of better agreement may 
be due to the fact that for Co solely experimental data for 
polycrystalline samples has been available. Remarkably, 
when the magnetization is changed from [001] into [110] 
direction, the sign change of a^^y in Co at the Fermi en- 
ergy level is compensated by the sign change of frjij,, and 

the net contribution stays the same. For Ni, the value 
calculated in bare GGA differs drastically from experi- 
ment, and even has the wrong sign. However, when the 
value of U is increased within GGA+C/, the calculated 
value approaches the experimental result in magnitude 
and sign. For U = 3.9 and J = 2.6, Q'/^+'j is about 87% 
of the smaller experimental value. This suggests that the 
main reason for the discrepancy between experiment and 
theory in Ni originates from an incorrect description of 
the Fermi surface within GGA. 



Overall, the values in Tab. (I) demonstrate that the 
intrinsic- and side-jump contribution play an important 
role in the ANE of ferromagnets. This finding is con- 
sistent with earlier studies in this field, which examined 
the behavior of the anomalous Nernst coefficient Qi as a 
function of the resistivity p^x- As far as the scattering- 
independent contributions are concerned, theory^^ pre- 
dicts a linear dependence of the form Qi/T cx pxx, which 
is also observed in experiment. '^^"^^ Our work substanti- 
ates these observations with quantitative analysis. 

The calculated temperature dependence of the ther- 
moelectric conductivity tensor in considered metals is 
depicted in Fig. (2). The importance of the side-jump 
contribution to the ANE is stressed by the fact that for 
the materials which are considered here a^^y is of the same 
order of magnitude or even larger than a^^y. In the zero 
temperature limit, one can apply the Sommerfeld expan- 
sion to the integral in Eq. (2) to obtain the standard Mott 
formula which relates the ANE to the energy derivative 
of the AHE:44 



3e 



da 



de 



(5) 



For the intrinsic Nernst effect, it is known that the above 
formula may be violated as T — >■ K.^ Indeed, we find 
that the energy derivative of the intrinsic contribution 
converges only very slowly with respect to the number 
of k-points in the Brillouin zone that are used for the 
evaluation of a^^y. The slow convergence is due to the 
sensitivity of the Berry curvature to the position of the 
Fermi energy, in particular when the latter approaches 
avoided band crossings or points of band degeneracy.^^ 
However, for the side-jump contribution, Mott's formula 
Eq. (5) holds, as can be seen by interpolating the curves 
in Fig. (2) to T = K. Apart from Ni, it yields a rather 
good estimate for the value of a^^y/T at room tempera- 
ture as well. 

In summary, we presented the ah initio calculations of 
the scattering independent contributions to the ANE in 
several ferromagnets. The theoretical values for the ther- 
moelectric conductivity tensor and the comparison of the 
calculated anomalous Nernst coefficient with experiment 
suggests that the ANE in elementary Fe, Co, Ni, in the 
ferromagnetic alloy FePt and presumably also in FePd 
is largely caused by the intrinsic and side-jump mecha- 
nisms. Discrepancies between theory and experiment in 
Ni are likely due to the imprecise description of correla- 
tion effects within bare GGA, which can be remedied by 
GGA-I-J7 calculations. 
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PACS numbers: 

We compare the theoretical value of the anomalous Nernst coefficient Qi = Pxx{oxy ~ Sa^y) /4:TtAI2 vi^ith 
experimental data. For this purpose, we need the values of the resistivity p^x, the thermoelectric conduc- 
tivity tensor Uxy, the Seebeck coefficient S, the conductivity axy and the magnetization M^. In linear 
transport theory, the transport coefficients axy = e£^°' and axy = C^^^ /T can be computed from the 
formula 

= -llde ^-^^^ ■ -xy ■ {e - pT . (1) 

where the expression for the Fermi distribution function is given by f{e,fi,T) = l/[exp((£ — ii)/kBT) + 1]. 
In Eq. (1), we insert the anomalous Hall conductivity (Jx^^ = '^^xy+'^xy for "'J^^j which has been computed 
with the same method as developed in Ref.^. In most experiments, the values of Qi, Pxx, S and are 
not measured simultaneously. As illustrated in Fig. (1), we therefore estimate pxx and S from linear 
interpolation whenever their values are not available, 

Pxx = ai-T + bi, S = a2-T + b2. (2) 

For AIz we take the saturation magnetization at room temperature. Since the Curie temperature Tc for 
Fe, Co and Ni is quite large, Tc 3> 300 K, we suppose that the temperature dependence of has no 
significant effect on the calculated values for Q^^'^^^ which are presented in Tab. (I). 
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FIG. 1: Interpolation of p^x and S for Fe, Co and Ni based on data from Ref.^ for polycrystalline samples. 
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TABLE I: Comparison with experimental data. Temperature T in K; magnetic field B in G; resistivity p^a: in 
/if2-cm; Seebeck coefficient 5 in /^V/K; anomalous Hall conductivity a^^y in l/(n-cm); anomalous thermoelectric 
conductivity a^^y in 10~'^A/(m-K); anomalous Nernst coefficient Q in 10~^^V/(K-G). Values that are labeled by 
an asterix are taken from interpolation of data from Ref.^"^. The sign of the values from Ref.*^ has been adjusted 
to be in accordance with the sign convention of Campbell*" . 
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Qi Expt. 


Fe 


291 


6290 


9.63* 


15.24* 


744 


0.505 


115 


— 1.254 


—0.43 


—0.24 


—0.67 


— 1.05''' 




293 




12.55 


14.29 


744 


0.509 


115 


— 1.254 


—0.52 


—0.30 


—0.82 


—0.82^ 




313 


23000 


11.65 


14.07* 


741 


0.544 


115 


— 1.258 


— 0.46 


—0.29 


—0.75 


—2.04** 




313 




13.73 


13.70 


741 


0.544 


115 


— 1.258 


—0.52 


—0.34 


—0.87 


—0.90'^ 




323 


23000 


12.21 


13.54* 


739 


0.559 


116 


— 1.260 


— 0.45 


—0.31 


—0.76 


—2.06* 




333 


18000 


12.49* 


13.01* 


738 


0.572 


116 


— 1.260 


—0.43 


—0.32 


—0.76 


—0.81^ 




333 


23000 


12.77 


13.01* 


738 


0.572 


116 


— 1.260 


-0.44 


—0.33 


-0.77 


—2.08** 




333 




14.91 


13.07 


738 


0.572 


116 


— 1.260 


-0.52 


—0.39 


—0.91 


—0.98^ 




343 


23000 


13.33 


12.48* 


737 


0.583 


116 


— 1.261 


—0.43 


—0.35 


—0.78 


—2.10* 


Co [onil 


297 


6290 


6.02* 


—30.09* 


484 


0.216 


222 


— 1.831 


0.51 


0.04 


0.55 


2.00^ 




297 


10570 


6.02* 


—30.09* 


484 


0.216 


222 


— 1.831 


0.51 


0.04 


0.55 


1.80^ 




298 


9000 


6.06* 


—30.18* 


484 


0.220 


222 


— 1.832 


0.52 


0.04 


0.56 


1.90" 




300 




6.13* 


—30.36* 


484 


0.226 


222 


— 1.833 


0.53 


0.04 


0.57 


1.80^^ 




320 


9000 


6.80* 


—32.16* 


485 


0.286 


223 


— 1.850 


0.63 


0.05 


0.68 


2.19" 




333 


18000 


7.24* 


—33.33* 


486 


0.320 


223 


— 1.861 


0.70 


0.05 


0.75 


2.00^ 




350 




7.81* 


-34.87* 


486 


0.359 


223 


-1.874 


0.80 


0.05 


0.85 


2.15^^ 


Co [1 1 nl 


297 


6290 


6.02* 


—30.09* 


125 


1.010 


—91 


1.657 


0.23 


0.07 


0.30 


2.00^ 




297 


10570 


6.02* 


—30.09* 


125 


1.010 


—91 


1.657 


0.23 


0.07 


0.30 


1.80^ 




298 


9000 


6.06* 


—30.18* 


125 


1.009 


—91 


1.644 


0.23 


0.07 


0.30 


1.90" 




300 




6.13* 


—30.36* 


126 


1.007 


—90 


1.618 


0.23 


0.07 


0.31 


1.80^^ 




320 


9000 


6.80* 


—32.16* 


128 


0.994 


—88 


1.359 


0.28 


0.06 


0.34 


2.19" 




333 


18000 


7.24* 


—33.33* 


130 


0.995 


—87 


1.197 


0.31 


0.04 


0.35 


2.00^ 




350 




7.81* 


—34.87* 


133 


1.006 


—85 


0.994 


0.36 


0.02 


0.38 




Ni Lf = 3.9 


291 


6290 


7.03* 


— 18.75* 


—842 


6.347 


444 


—0.379 


0.33 


0.90 


1.23 


3.555 




291 


10620 


7.03* 


— 18.75* 


—842 


6.347 


444 


—0.379 


0.33 


0.90 


1.23 


1.30^ 




300 




7.42* 


-19.19* 


-843 


6.330 


443 


-0.355 


0.37 


0.97 


1.34 


5.25" 




311 


9000 


7.90* 


-19.73* 


-844 


6.307 


443 


-0.327 


0.41 


1.08 


1.49 


2.59" 




313 


6670 


12.15 


-19.83* 


-844 


6.302 


442 


-0.322 


0.64 


1.66 


2.30 


5.93* 




323 


6670 


12.63 


-20.32* 


-845 


6.276 


442 


-0.298 


0.69 


1.79 


2.48 


6.62* 




330 


9000 


8.72* 


-20.66* 


-846 


6.257 


441 


-0.282 


0.49 


1.26 


1.75 


3.04" 




333 


6670 


13.04 


-20.81* 


-846 


6.248 


441 


-0.275 


0.74 


1.90 


2.64 


7.31* 




343 


6670 


13.46 


-21.30* 


-847 


6.218 


440 


-0.253 


0.78 


2.02 


2.80 


8.00* 




350 


9000 


9.59* 


-21.64* 


-847 


6.195 


439 


-0.237 


0.57 


1.47 


2.03 


3.66" 




350 




9.59* 


-21.64* 


-847 


6.195 


439 


-0.237 


0.57 


1.47 


2.03 


7.25" 


FePt 


293 


52566 


35.90 


-12.40 


828 


2.007 


132 


0.957 


4.61 


1.27 


5.88 


5.60^^ 



Fe 4nM^ = 22100 G* 

Co 47rM, = 17857 G'^ 

Ni 47rM, = 5670 G* 

FePt 4nM, = 12566 G^^ 



